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SMALL-TIME  BEHAVIOR  OF  UNSTEADY  CAVITY  FLOWS 


by 

D.  P.  Wang  and  T.  Yao-tau  Wu 
California  Institute  of  Technology 

A  perturbation  theory  ie  applied  to  inveetigate  the  email-time  behavior  of 
unsteady  cavity  flows  in  which  the  time -dependent  part  of  the  flow  may  be 
taken  as  a  snoall-time  expansion  superimposed  on  an  established  steady 
cavity  flow  of  an  ideal  fluid.  One  purpose  of  this  paper  is  to  study  the 
effect  of  the  initial  cavity  sise  on  the  resulting  flow  due  to  a  given  dis¬ 
turbance.  Various  existing  steady  cavity-flow  models  have  been  employed 
for  this  purpose  to  evaluate  the  initial  reaction  of  a  cavitated  body  in  an 
unsteady  motion.  Furthermore,  a  physical  model  is  proposed  here  to 
give  a  proper  representation  of  the  mechanism  by  which  the  cavity  volume 
may  be  changed  with  time;  the  initial  hydrodynamic  force  resulting  from 
such  change  is  calcxilated  based  on  this  model. 


1.  Introduction 

It  has  been  noted  that  the  free  surface  flow  problems  become  very 
difficult  to  treat  exactly  when  the  flow  is  time-dependent.  Some  of  the 
principal  difficulties  involved  in  this  class  of  problems  can  be  envisaged 
as  follows.  First  of  all,  the  problem  is  generally  nonlinear  because  the 
boundary  conditions  on  the  unsteady  free  surface  are  nonlinear  and  be¬ 
cause  the  location  of  the  free  boundary  is  not  known  i,  priori.  Secondly, 
unlike  the  case  of  steady  flows  with  free  boundary,  the  surface  of  constant 
pressure  is  no  longer  a  surface  of  constant  speed.  Consequently,  tiie 
powerful  hodograph  method  widely  used  in  tiie  steady  free  boundary 
probleihs  loses  nwch  of  its  potential  and  usefulness  for  possible 
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applications  in  this  case.  Probably  for  this  reason,  only  a  relatively  small 
number  of  unsteady  free  surface  problems  (of  racher  special  cases)  have 
been  treated.  Some  comprehensive  survey  and  review  of  the  recent 
literatures  are  available  elsewhere  (see  Cilbarg[l],  Birkoff  and 
Zarantonello[  2]  and[3]  ),  and  hence  will  not  be  recited  here. 

It  is  worth  noting  that  there  is  an  essential  difference  between  the 
unsteady  flows  with  and  without  a  free  surface.  In  the  determination  of 
the  velocity  field  of  an  unsteady  flow  without  a  free  surface  (which  includes 
the  vortex  sheet  shed  from  an  oscillating  wing),  the  time  appears  only  as 
a  parameter.  The  kinematics  of  the  unsteady  motion  will  therefore  not 
differ  basically  from  its  corresponding  steady  flow.  On  the  other  hand, 
for  an  unsteady  flow  with  a  free  surface,  the  time  variation  affects  the 
flow  explicitly  through  the  boundary  conditions  on  the  free  surface.  There¬ 
fore,  the  flow  will  depend  on  all  its  previous  time  history,  which  further 
complicates  the  detailed  analysis  and  calculation.  However,  when  an 
established  steady  cavitating  flow  (regarded  as  the  basic  flow)  is  given  a 
sudden  acceleration,  the  problem  of  the  small-time  behavior  of  the  flow 
is  relatively  simple  because  there  will  be  only  a  short  history  of  auiy  time- 
varying  disturbance  to  be  considered  in  this  period.  It  is  therefore  clear 
that  some  difficult  problems  in  unsteady  cavity  flows  can  be  best  handled 
in  the  small-time  limit  since,  for  one  reason,  the  auialysis  is  relatively 
straightforward.  We  hope  that  such  an  investigation  may  cast  light  on 
some  of  the  basic  features  of  the  unsteady  cavity  flow,  especially  the  be¬ 
havior  of  the  large  frequency  limit  in  the  case  of  oscillating  cavity  flows. 

The  impact  problem  of  the  Kirchhoff-Helmholtz  flow  with  a  flat 
plate  normal  to  the  stream  has  been  treated  by  Gurevich[4].  Attempt  is 
made  here  to  generalize  the  consideration  to  include  the  effects  of  (1)  the 
cavitation  number  (the  cavity  size  being  finite),  (2)  the  angle  of  attack, 
and  (3)  the  time  variation  of  the  net  cavity  volume. 

The  existing  mathematical  theories  for  steady  cavity  flows  are 
seen  to  be  based  on  some  proposed  physical  models.  It  is  well  known  that 
in  a  positive  range  of  the  cavitation  number  up  to  moderate  values  of  order 
unity  the  agreement  between  the  existing  models  for  steady  flows  may  be 
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considered  very  clo8e[5].  No  comparison,  however,  has  ever  been  made 
(to  the  authors'  knowledge)  between  these  existing  models  for  the  case  of 
unsteady  cavity  flows.  In  this  paper,  the  re-entrant  jet  model, 
Riabouchinsky's  model,  and  a  modified  wake  model  are  applied  to  in¬ 
vestigate  the  small-time  behavior  of  the  unsteady  cavity  flow,  including 
the  effects  of  the  cavitation  number  and  the  incidence  angle  of  the  flow. 

The  problem  related  to  the  possibility  for  the  volume  of  a  finite  or 
infinite  cavity  to  vary  with  time  deserves  a  particular  note.  For  cavity 
flows  of  an  incompressible  liquid  surrounding  a  vapor  cavity,  it  is 
obvious  that,  when  the  cavity  volume  changes,  the  conservation  of  mass 
and  conservation  of  volume  of  the  entire  flow  becomes  incompatible  be¬ 
cause  of  the  difference  in  the  densities  of  these  two  phases.  Consequently, 
any  variation  of  the  cavity  volume  must  come  from  a  source  distribution 
in  the  flow  with  its  net  strength  depending  on  the  time -rate  of  change  of 
the  cavity  volume.  It  has  been  proposed  by  one  of  the  authors [ 6]that  as  a 
physical  model,  the  change  in  the  cavity  volume  can  be  affected  by  a 
source  or  sink  located  at  the  point  of  infinity.  A  direct  consequence  of 
this  source  with  a  time-varying  strength  in  a  two-dimensional  flow  of 
infinite  extent  is  that  it  generates  a  pressure  field  which  is  logarithmical¬ 
ly  singular  at  infinite  distances.  To  create  such  a  flow  it  requires  an 
infinite  amount  of  energy.  Probably  for  this  reason  the  treatment  of  the 
unsteady  partial  cavity  flow  by  Timman[7],  the  work  of  Geurst[8],  of 
Parkin[9],  and  of  Wu[6],  have  been  carried  out  under  the  assumption  of 
fixed  cavity  volume  and  no  source  at  infinity.  In  reality,  however,  the 
flow  is  usually  finite  in  extent  and  never  two-dimensional  in  the  large;  the 
pressure  singularity  at  infinity,  being  outside  the  flow  region,  may  there¬ 
fore  be  regarded  as  a  simplifying  idealization,  much  the  same  as  the 
representation  of  a  solid  body  by  a  distribution  of  singularities.  This 
physical  model  for  permitting  the  cavity  volume  to  change  has  been 
applied  to  various  cases  together  with  other  cavity  flow  models.  It  is 
expected  that  this  variable -cavity- volume  model  will  give  a  good  representa¬ 
tion  of  the  flow  field  near  the  body- cavity  system.  The  validity  of  this 
model  of  course  will  have  to  be  verified,  directly  or  indirectly,  by  future 
experimental  observations. 
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2.  General  Formulation 

To  fix  idea  ,  we  suppose  that  for  the  time  t  <  0,  a  steady  plane 
flow  of  an  incompressible,  inviscis  fluid  past  a  solid  body  with  a  cavity 
formation  has  been  established,  the  solution  of  which  is  assumed  to  be 
given,  or  can  be  determined  with  the  aid  of  some  cavity  flow  models. 
Suppose  now  the  solid  body  to  which  the  cavity  is  attached  is  given  for 
t  >  0  a  sudden  accelerated  motion;  the  problem  is  to  evaluate  the  small¬ 
time  behavior  of  the  resultant  unsteady  cavity  flow. 

In  general,  the  motion  of  the  rigid  boundary  may  consist  of  a 
translation  and  a  rotation.  Let  (x^,  y^)  be  a  point  on  the  rigid  surface 
S^(x^,  y^)  =  0  in  the  basic  steady  flow,  and  let  it  be  displaced  in  time  t 
to  the  position  (x,  y)  with  translational  velocity  (V^{t),  V^(t)  )  and 
angular  velocity  u(t).  In  terms  of  the  complex  variable  z  =  x  iy  and 
V(t)  =  +  i  V^,  the  motion  of  z  may  be  written 

dz/dt  =  V(t)  +  iw(t)  z  .  (1) 

We  shall  assume  that,  for  small  positive  t,  V  and  u  may  be  expanded 
in  power  series  of  t,  starting  with  the  linear  term  so  that 


dz/dt  =  c  t  +  c  t*  +  i  (w  t  +  w  t*)  z  +  0  (t^ )  (2) 

12  12 


where  c  =a  +ib,a,b  and  w  being  real  constants.  It  then 
n  n  n  n  n  n  " 


follows  that  for  small  t, 


+  i  (c  +iwz^)t*+  4  (c  +iwz  )t*+  0(t^)  (3) 


Z  =  Z_1--5\C  TIW  V  T  XV,  T  *-  -  , 

2  I  I  o'  3  j  j  o' 


The  displaced  surface  will  be  denoted  by  S(x,  y,  t)  =  0.  In  fact,  we  have 


So(xQ{x,y,t)  ,  )  =  S(x,y,t)  , 


(4) 
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regarding  (1)  as  to  provide  the  canonical  transformation  =  x^(x,y,t), 

Vo* 

From  the  nature  of  the  body  motion  it  also  follows  that  the  com¬ 
plex  velocity  potential  of  the  flow 


f(2,t)  =  i;(>(x,y,t)  +  i4(  (x,y,t)  ,  (5) 


will  assume  for  small  t  the  expansion 

f(2.t)  =^2)  +  tf^{*)  +  4  +  ••••  .  (6) 

Where  f„(*)  =  <»„(x,y)  +  i«|»Jx.y).  n  =  0,1,2,  and  r{2)  is  the 

complex  velocity  potential  of  the  basic  flow.  The  function  (x,y)  may 
be  called  the  initial  acceleration  potential.  While  4;^  =  const,  gives  the 
streamlines  of  the  basic  flow,  the  harmonic  functions  for  ni^l, 
being  the  complex  conjugate  of  <p^,  are  introduced  solely  to  make  f^(z) 
analytic  functions  of  z.  The  velocity  components  are,  as  usual 

u(x,y,t)  =  d<p/dx  ,  v(x,y,t)  =  dip/ &y  .  (7a) 

We  shall  introduce  the  complex  velocity  w  =  u  -  iv,  and  w^  =  df^/dz, 
then  from  (6) 

w(z,t)  =  w^(z)  +  tw^(z)  +  ^  t*w^(z)  +  . .  .  .  (7b) 

Similarly  the  pressure  p(x,y,t)  may  be  assumed  to  possess  the 
expansion 

p(x,  y,  t)  =  p^(x,  y)  +  p^(x,  y)  +  t  p^(x,  y)  +  t*p^(x,  y)  +  . . .  .  ,  (g) 

where  p^  denotes  the  pressure  field  of  the  basic  flow,  p  the  impulsive 
pressure  due  to  the  sudden  acceleration.  Then  from  the  Bernoulli  equation, 

p/p  +  dip/dt  +  ^  (V«)*  =  p^/p  +  .^  U*  (9) 
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where  p^and  U  are  the  pressure  and  speed  of  the  basic  flow  at  infinity, 
we  obtain,  by  equating  the  coefficients  of  same  powers  of  t,  the  following 
relations: 

Pq/p  -Z  =  Poo^P  +  4 

Pj/p  =  ■ 

■p  ! p  =  ~  <p  -  9  <p  )  ,  and  so  on. 

The  boundary  conditions  of  the  problem  are  as  follows: 

(i)  At  the  solid  surface,  the  normal  component  of  the  flow  velocity 
relative  to  the  moving  boundary  must  vanish.  An  alternative  way  of 
stating  this  condition  is  that  the  fluid  particles  originally  on  S(x,y,t)  =  0, 
at  small  time  interval  apart,  will  remain  on  it.  That  is 

DS  _  SS  .  8^  8S  8^  8S  _  0  S*  v  t^  -  0 

Tff-  -  +  +  5(x,y,t)  -  U  , 


which  becomes,  upon  using  (4), 


®So 

'b^  ^  8^  ®*o  ^ 

BS 

.  0 

a  By  a  By  By  ' 

B9  '0  B9  '0  ^0 

~5r 

0 

•55^  55r  -57  -Sy-  W- 

K 

’0 

w  ^  ^  "St” 

>  ¥ 

j 

(H) 


on  S(x,y,t)  =  0.  Here  the  functions  x^(x,y,t)  and  y^(x,y,t)  can  be 
written  down  immediately  from  (3)  by  interchzmging  z  and  z^,  and  by 
changing  the  signs  of  a^,  b^  and  as  can  readily  be  seen  from  (3), 

Equation  (11)  is  in  a  form  convenient  for  manipulation  since  8S^/8x^  and 
8S^/8y^  already  correspond  to  the  components  of  the  normal  to  the  initial 
surface  ~  Substituting  (3)  and  (6)  in  (11),  expanding  the 

various  quantities  about  (x^,y^)  and  t  =  0,  and  equating  the  coefficients 
of  different  powers  of  t,  we  obtain  the  conditions  that  on  S^(x^,y^)  =  0, 

d9  /Bn  =  0  , 
o  o 


(12a) 
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I  /  0n  =  n  (a-u  y  )  +  n  (b  +  w  x  )  , 

1  o  o  I  1  o  o,  »  »  ® 

1  2 


(12b) 


(8®  8^ 
n  TE—  -  TT- 

o  By  o  ox 
.  1  2 


J.  = 

1 

o  1 


o 

w: 


-  (a  -w  y  ) 
1  1^0 


02- 


O  2 


I  8x 


n  -B — -(b  +<o  X  ) 
o  oxoy  '  1  I  o' 
2  o 
■*  ^ 


8*. 


2-  1 


S  .  'O 


aV 


”o  8x  8y  ft,  2 

1  o  'o  2  8y^ 

(12c) 


where  n  =  (n  ,  n  )  is  the  unit  outward  normal  to  the  surface 
o  o  o 


(ii)  There  are  two  boundary  conditions  at  the  free  surface  of  the 
cavity.  Suppose  the  free  surface  may  be  expressed  as 


F(x,y,t)  =  y  -  h(x,t)  =  0  , 

then  the  kinematic  condition  that  the  fluid  particles  on  the  free  surface 
will  remain  on  it  requires 

V  =  8h/  8t  +  u  8h/  8x  on  y  =  h(x,  t)  .  (13) 

We  assume  that  for  small  t,  h(x,  t)  may  be  expanded  as 

h(x,t)  =  h  (x)  +  t  h  (x)  +  4  t*h  (x)+o(t^)  ,  (14) 

u  1^2 


where  y  =  h^(x)  denotes  the  cavity  boundary  of  the  basic  flow  on  which 

dh^/dx  =  v^(x,h^)/u^(x,h^)  .  (15) 

Substituting  (7)  and  (14)  in  (1 3),  and  expanding  u^^  and  v^  on  y=h(x,t) 
about  y  =  h^(x)  and  t  =  0,  we  obtain 

\  =  V  <*^0/*^  =  ®  ’ 

h^  =  v^(x,ho)  -  Ui(x,ho)<iho/dx  . 


(16) 
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Since  s  0,  the  free  eurface  will  not  be  displaced  in  the  time  of  order 
t,  as  should  be  expected. 

In  the  basic  flow,  the  cavity  boundary  is  a  surface  of  constant 
pressure,  and  hence  also  one  of  constant  velocity,  say 


Po'Pc  ’  '“Ic  ’  y  •  (17) 

so  that  by  the  Bernoulli  equation  (10), 

q^=U(l+a)*  (18a) 

where  a  is  the  cavitation  number  of  the  basic  flow,  defined  as 

o  =(P«-Pc)  /  -Z  PU*  • 

We  shall  assume  that  the  pressure  in  the  cavity  of  the  unsteady  flow  will 

be  maintained  at  the  same  constant  value  p  ,  that  is, 

^  c 

p(x,y,t)=p^  on  y=h(x,t)  .  (19) 

By  expanding  the  left  side  of  this  dynamic  condition  about  y  s  h^(x)  and 
t  =  0,  using  (8),  (14)  and  (16),  the  following  conditions  result 

p  =p  =0  on  y=h  (x)  .  (20) 

12 

Hence  from  (10), 


=  0  on  y  =  h^(x)  , 

9  ^  =  ony=h^(x)  . 


(21a) 


Upon  differentiating  (21a)  along  y  =  h^(x)  and  using  (15),  it  is  readily 
seen  that  (V^^)  •  (V ^^ )  =0,  and  hence 


=  0  on  y  =  h^(x)  . 

From  (21a)  and  (21b)  it  therefore  follows  that 

(V<»^)‘(Vf^)  s  ('^9^)'{y<P^)  =  0  on  y  =  h^(x) 


(21b) 


(22) 
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which  asserts  that  the  perturbation  velocity  is  normal  to  the  original  cavity 
boundary  up  to  the  time  of  order  t^. 

(iii)  At  the  point  of  infinity  we  require  the  perturbation  velocity  to 

vanish, 

I  V^n  I  ““*■  0  I  2  I  — *•  00,  for  n  =  1 , 2,  •  •  •  .  (23) 

If,  in  addition  to  the  sudden  acceleration  of  the  solid  body  and  the  assump¬ 
tion  of  the  constancy  of  the  cavity  pressure,  a  certain  amount  of  fluid  is 
removed  at  infinity  with  a  source  strength  +  •  •  •  so  that  the 

cavity  volume  may  be  changed  arbitrarily,  then,  aside  from  condition  (23), 
we  must  impose  additional  conditions  at  infinity  as 


Im 


n  =  1, 2 ,  *  •  •  , 


(24) 


where  r  is  a  contour  around  the  point  of  infinity.  However,  we  impose 
no  boundedness  condition  on  p  at  infinity  since,  if  such  can  be 
arbitrarily  chosen,  p^(n  >  1)  will  be  logarithmically  singular  at  infinity. 
It  will  be  seen  later  that,  when  the  cavity  is  taken  to  be  infinitely  long 
(Helmholtz  flow,  corresponding  to  the  cavitation  number  v  =  0),  the 
solution  exists  only  when  =  0.  Consequently  the  effect  of  change  in 
cavity  volume  can  be  sought  only  in  the  case  of  finite  cavities.  However, 
the  limit  of  the  hydrodynamic  forces  in  such  cases  as  the  cavity  becomes 
infinitely  long,  with  held  fixed,  is  seen  to  exist. 

On  the  other  hand,  from  Kelvin's  theorem  on  the  conservation  of 
circulation,  the  circulation  around  the  point  of  infinity  cannot  be  changed 
in  the  unsteady  motion  for  t  <  so.  By  combining  this  condition  with  (24), 
we  may  write 


n  =  1 , 2,  3  •  •  •  . 


(25) 


This  completes  our  formuleition  of  the  small-time  perturbation  theory. 
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3.  Solution  of  the  Perturbed  Flow. 

For  the  moment  we  assume  that  the  solution  of  the  basic  steady 
flow  is  given.  We  note  that  in  the  basic  flow  the  entire  boundary  of  the 
body- cavity  system  belongs  to  a  streamline,  the  form  of  which  is  known. 
Therefore,  by  the  general  theory  of  conformal  mapping,  it  is  always 
possible  to  find  an  analytic  function 

4  (z)  =  S(x,y)  +  i  ii(x,y)  (26) 

such  that  the  entire  flow  region  in  the  z-plane  is  mapped  by  (26)  into  the 
upper  half  ^-plane  with  the  entire  boundary  of  the  body-cavity  system 
mapped  onto  the  entire  |-axis.  For  simplicity,  we  shall  make  the  part 
|4l  1  of  the  real  4~^xi8  correspond  to  the  wetted  solid  surface  and  the 

part  III  >1  to  the  cavity  boundary. 

After  the  transformation  to  the  I'ple'Oe  the  boundary  value  problem 

formulated  in  the  last  section  can  be  stated  as  a  Hilbert  problem,  as  will 

be  shown  below,  the  solution  of  which  is  readily  attained.  It  is  noted  from 

the  last  section  that  the  problems  of  different  orders  in  the  perturbation 

are  expressed  in  a  similar  form.  That  is,  the  normal  velocities  d«  /Bn 

no 

are  given  at  the  initial  solid  surface  and  the  potentials  are  prescribed 

at  the  unperturbed  cavity  boundary.  It  is  therefore  sufficient  to  treat  one 

as  the  typical  case.  To  save  writing  we  shall  denote  f  =  v  +  iU*  bv 

n  ’^n  ^n  ^ 

F  =  *  +  i 


It  is  convenient  to  introduce  the  analytic  function 
G(4)  =  dF/d4  =  (dF/dz)  (dz/d4) 


(27) 


which  is  defined  for  Im  4  >  0.  Then  on  the  solid  surface,  n  =  0+,  |||  <  1, 


Im  G  = 


8« 


8ft 

•RT 


dz 

HI 


(6) 


say. 


(28) 


which  is  known,  by  (12)  and  (26).  Furthermore,  on  the  cavity  boundary, 

T1  =  0  +  ,  III  >  1, 
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Ri  G  =  8»/8e  =  g  (fi)  say. 
2 


(29) 


which  is  also  known,  by  (21).  In  particular,  g  =  0  for 

2 

Finally,  (25)  becomes 


Ir  JE  ir,  = 


and  <p^ . 


(30) 


where  F'  is  the  image  of  F  in  the  ^plane,  the  subscript  n  of  Q  is 
dropped  to  show  a  typical  case. 

It  is  possible  to  transform  the  above  boundary  value  problems  of 
a  mixed  type  into  a  Hilbert  problem  by  extending  the  unknown  function 
G(4)  to  a  sectionally  analytic  function,  defined  on  the  whole  (-plane 
(excluding  the  real  (-axis  if  necessary).  Since  ^  and  higher  terms  will 
not  be  treated  explicitly  here,  we  shall  demonstrate  the  method  by  taking  g^ 
of  (29)  to  be  zero;  the  general  case  of  g^  /  0  can  be  carried  out 
similarly.  First,  the  function  G(()  may  be  continued  into  the  lower 
half  (-plane  by 


G(X)  =  -  GT?:)  .  (31) 

For  the  case  g  : 

2 

the  interval  |  (  | 
limiting  value  of 
follows  that 

G^  +  G’  =  2i  Im  G^  =  2i  g^  (()  for  |  (  |  <  1  , 

G’*' -  G"  =  2  Ri  G’*’ =  0  for  1(1  >1  .  (32) 


9,  G(()  is  the  analytical  continuation  of  G(()  through 
>  1.  In  the  following,  G-(()  will  be  used  to  denote  the 
G(()  as  ►  +  0.  From  (28),  (29)  and  (31)  it  then 


The  above  Hilbert  problem  is  well-known  [lO],  its  general  solution  cam 
be  written 


G(()  = 


(t*-l) 


-i 


i; 


(i-e‘) 

“TT 


g,(e) 


d(  + 


iP(C) 


(;*-i) 


\ 


for 


n  >  0 

(33) 
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where  the  function  (C*  -  1)  is  defined  on  the  entire  (-plane  with  the 
branch  80  chosen  that  ((*  -  1)^“*(  as  !(!“♦«,  and  P(()  is  an 
arbitrary  Laurent's  series  with  real  coefficients.  The  last  term  in  (33)  is 
the  general  solution  to  the  homogeneous  problem  (with  g  =  0  also).  The 
real  coefficients  of  P((),  and  hence  G((),  can  be  determined  uniquely, 
when  (30)  is  satisfied  and  the  condition  that  the  pressure  is  integrable 
over  the  rigid  boundary  is  observed. 

In  the  following  the  above  perturbation  theory  will  be  carried  out 
for  several  basic  steady  cavity  flows. 


4.  Inclined  Lamina  in  Kirchhoff  Flow 

As  a  simple  example  we  consider  the  basic  flow  to  be  that  past  a 
flat  plate  held  at  an  &ngle  a,  with  a  cavity  formation  of  infinite  length. 
The  solution  of  this  problem  is  known  [11],  which  we  simply  cite  below 
for  the  subsequent  use.  The  coordinate  system  in  the  z-plane  and  its 
conformal  mapping  planes  are  shown  in  figure  1.  For  simplicity  the 
free  stream  velocity  U  and  the  plate  length  i  ,  are  normalized  to  unity. 
The  solution  w^  =  w^(z^)  can  be  written  parametrically  as 


w 


=  (/r7p  -!)/(  , 

=  -Kj^  (  JTTT*  +  1)  ((  + sec  a)-*  dt 


K 

7 


sin  a  ((  +  sec  a) 


((  +  sec  «)* 

% 

where 

K  =  2  tan’  or  sin  al{4  +  ir  sin  a)  . 


TK 


(34a) 


tan  -  (l+sec  or) 
(34b) 

(34c) 


The  entire  flow  region  is  mapped  by  (34)  into  the  upper  half  (-plane  with  the 

corresponding  boundary  as  prescribed  in  the  last  section 

(i.e.  on  T)  =  0,  |(  |<  1  corresponds  to  the  plate  and  |  (  |  >1  to  the 
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cavity  surface).  The  coefficient  of  the  normal  force  of  the  basic 
flow  is 

C..,  =  N  /(  i  pU*  i)  =  2ir  sina/(4  +  ir  sina)  .  (35) 

O  ^ 

O 

(a)  The  first  order  solution 

The  unit  normal  to  the  plate  is  now  n^=(0,  -1).  We  suppose  the 
rotational  motion  is  referred  to  the  leading  edge  of  the  plate.  Then  the 
boundary  conditions  (12b),  (21a)  and  (23)  become 

/  0y  =  (b  +w  X  )  on  =  0  ,  |  ^  |<  1  ,( 36a) 

=  0  and  hence  b<pj  8g=0  onii  =  0,  16|>1.  (36b) 

(V^  I - -  0  as  t - -  -  seca  ,  (36c) 

s 

It  is  noted  that  a^,  the  first  term  of  the  x-component  of  acceleration, 

drops  out  from  (36a), implying  that  the  acceleration  of  the  plate  parallel  to 

itself  has  no  effect  on  the  flow  up  to  time  of  order  t.  Furthermore, (25) 

cannot  be  satisfied  unless  all  the  Q  =0.  This  can  be  seen  as  follows. 

n 

If  f  0,  then  (25)  implies  that  w^  =  Q^lZvz  +  o()z|"*),  and  hence 
~(Q^/2ir)  log  |z|,  as  |z| — ►«.  It  follows  that  p^  will  be  logarith¬ 
mically  singular  at  z  =  oo  ,  which  contradicts  the  conditions  of  p^,  such 
as  (20),  on  the  undisturbed  free  boundary  y  =  h^(x)  which  extends  to  in¬ 
finity.  This  indicates  that  the  Kirchhoff-Helmholtz  model  with  an  infinite 
cavity  is  not  a  realistic  model  for  the  consideration  of  change  of  cavity 
volume.  The  problem  when  ^  0  will  be  considered  later  when  other 
finite-cavity  models  are  adopted. 

By  making  use  of  (36a),  (36b)  and  (34b),  (28)  and  (29)  become 
ImG=^  I  ^  I  =K(b^+(ox^)(/iTp+l)(e  +  8eco)“*  .  |g|<  1 


Rf  G  =  0  , 


le  l>  1 
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where  x^(4)  i*  given  by  (34b),  and  K  ia  given  by  (34c).  Finally,  by 
(33),  we  obtain 


df 

_ 1  _ 

ST 


j.  (pi  (b +«  X  )(i-e*+JTl5^) 

— _!_-L.2 - L.^ -  dt  +  )  C  (l+8ec«)" 

<S-t)(l+*ec«)  J 

(37) 


where  are  real  coefficients.  The  Laurentz  series  in  the  last  term  is 
expanded  about  4  =  -  seca(or  z  =  oo)  for  the  convenience  of  application 
of  the  boundary  conditions  at  z  =  cc.  The  first  term  in  (37)  behaves  like 
4  as  (4l  00  and  is  regular  everywhere  except  at  4  =  +  1.  Now 

by  using  (34b),  ~ 

-  ‘it  -  _  (4-t-geca)^ 

"  '  K(jr:4'+i)  ^  • 


w 


d4  Ha 


(38) 


From  ( 37)  and  ( 38)  we  readily  see  that  ®  0  except  for  n  =  -1,-2 


in  order  that  w  — ►  0  as  z 

t 


00  (or  4 


regular  as  4  — *■  *o- 
(30)  requires  that 


-  sec  a)  and  w^  is 


Furthermore,  with  =  0,  as  explained  above. 


zw  0 
1 


as 


(39) 


From  (34b),  we  find  that  as  (z  ( — , 


^  K  1  +  i  tan  a 

X  _  - 

(4  +  sec  a)* 


Making  use  of  this  result,  we  readily  deduce  from  (37)  and  (38)  that  (39) 
is  satisfied  if,  and  only  if,  C  =  C  =0  also.  Therefore  all  the  co- 

-I  -2 

efficients  vanish,  thereby  the  first  order  solution  given  by  (37)  and 
(38)  is  uniquely  determined. 

It  may  be  noted  that  w^  has  a  singularity  at  the  edges  of  the 
plate,  physically  corresponding  to  a  narrow  spray  sheet.  For,  with 
4  =  +  (1+e),  Ic  I  <  <  1,  we  deduce  from  (37)  and  (38)  that  as  4  +  1. 


4 
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(2«)* 


+  i 


(b^+ox^Kn-ji-eMj 


(S  +  sec  Of)* 


de  [i+o(|€  1^)] 


From  the  behavior  of  on  the  free  surface  (c  positive  small)  it  may 
be  seen  that  the  free  surface  starts  to  move  into  the  cavity  when  the 
plate  accelerates  into  the  fluid  (e.g,  ,  with  b  <  0,  =  0),  and  vice 

versa.  Furthermore,  it  is  noted  that  w  is  of  order  for  large 

values  of  (  z  |.  The  net  effect  can  be  seen  to  be  such  that  there  will  be  no 
net  change  in  the  cavity  volume. 

Since  the  spray  sheets  do  not  produce  any  singular  force  (unlike  the 
leading  edge  suction  on  a  thin  airfoil),  the  normal  force  acting  on  the 
plate  can  be  obtained  by  integrating  the  pressure  along  the  plate  so  that 
for  small  t,  in  view  of  (3), 


=  N  +  N  +  t  N  +  0(t*)  . 

O  1  2  ' 


(40a) 


The  first  term  is  given  by  (35).  Now,  from  (10) 

^  Pi  =  -p  £  0)^0  =  -pj  j 


Xo(e) 


»»;(g.o) 


dg  . 


(40b) 


By  substituting  the  real  part  of  (37)  for  9^  / 8g  in  the  above  integral,  the 
normal  force  coefficient  may  be  expressed  as 


(41a) 


where 
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the  above  integrals  being  interpreted  by  their  Cauchy  principal  values. 

In  (41a),  b^  and  are  expressed  in  the  physical  units,  and  the  plate 
length  i  and  free  stream  velocity  U  are  restored  for  completeness. 

The  integrals  in  (41b)  and  (41c)  cannot  be  expressed  in  terms  of  elementary 
functions  of  or:  they  are  integrated  numerically  and  the  results  are 
plotted  versus  a  in  figure  5.  In  particular,  we  have 

(ir/2)  =  0.8448  ,  (42) 

1 

which  is  the  special  case  treated  by  Gurevich  [4].  The  quantity 

1 

represents  the  jump  in  the  normal  force  coefficient  due  to  the  accelera¬ 
tion. 

(b)  The  second  order  solution  for  osir/Z,  u  =b  =c>)  =0  , 
_ lit. 

To  facilitate  investigation  of  the  behavior  of  the  second  order 
solution,  let  us  choose  the  special  case:  a  -  r!  Z  (the  flat  plate  being  held 
normal  to  the  stream)  and  b^,  all  vanish.  Then  the  boundary  con¬ 
ditions  (12c),  (21  ),  and  (23)  become 


8^^/8yo  = -bi8*W9yo^  , 

on  T)  =  0  ,  1  <  1  . 

(43a) 

ISI 

II 

O 

on  q  =  0  ,  i  £  1  >  1  , 

(43b) 

|V,  I  0 

as  |z  1  — ►  so  . 

(43c) 

Since  the  component  of  the  acceleration  parallel  to  the  plate  does  not 
appear  in  the  above  conditions,  the  flow  will  not  be  affected  by  it  up  to  the 
second  order  terms.  Now,  in  the  limit  as  a — *■  ir/ 2,  the  zeroth  order 
solution  becomes 

Zq  =  [  2(l-t)  (1-t*)^  +  cos'*!il/(4+tr) 

and  w^  is  still  given  by  (34a).  Consequently,  from  (28)  and  (29) 


(44) 
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8u 


Im  G  =  b 


|8x 


'SI 

Ri  G  =  0 
Hence,  by  (33), 


b^TMd-e*)'^- 1] .  Ti  =  0  ,  lei  <  1  . 

Ti  =  0  ,  I e l>i  . 


^  =  /  r  ‘  d{ 

^  irfi^  *^-1  (e-oe* 


(45) 


in  which  the  complementary  solution  is  zero,  as  can  be  shown  by  the 
same  argument  given  previously  for  <p  .  Carrying  out  the  integration,  we 
find 


w 

2 


df 

_ _ 2  _ 

-HI”  - 


(1-1*)“^ 

i(4-Hr) 

4t* 


{21  + log  -  iiKi-;")^}, 

i  df 

[i-(i-;vi^  . 


(46a) 

(46b) 


It  is  readily  verified  that  w^  also  behaves  like  for  large  |z  |  , 


The  second  order  normal  force,  by  (40a)  and  (10),  can  be  obtained 


from 


(  ^  I  ^  o 

*”  =  0  J.j  I  »/  “o".),,.  0  W 


=  P 


8^  dx 

ST 


de 


After  substituting  the  various  terms  and  evaluating  the  resulting  integral, 
we  find  the  simple  result: 

N  =  0  .  (47) 

2 

For  this  special  case  we  therefore  have  the  normal  force  coefficient 
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Cj,(t)  =  Cj^  +  Cj,  +  Cj^t+  0[(Ut/i)*) 

O  »  2 

=  0.8798  +  0.8448  (ib^/U*  )  +  0  (Ut/i)*  .  (48) 

Thus  for  constant  acceleration,  has,  aside  from  the  stepwise  change, 

the  following  behavior 

However,  N  may  not  vanish  when  b  ,  u  ,  and  u  are  different  from 
2  2  12 

zero. 

The  quantity  m^  =  ^ may  be  called  the  initial  induced  mass 
of  this  cavity  flow ,  then  by  (48) 

m  =  -^  pU*iCj^  b'  =  0.4224  pf*  .  (49) 

i  1 

If  the  flat  plate  had  undergone  an  acceleration  b^  normal  to  the  flow 
without  wake  formation  (a  postulated  Dirichlet  flow),  then  the  induced 
mass  would  be 


m^  •*  irp(//2)*  . 


(50a) 


Thus 


m  Im  =  0.5377  . 

1  o 


(50b) 


This  ratio  is  less  than  unity,  as  should  be  expected,  since  the  cavitated 
side  of  the  plate,  being  exposed  to  constant  pressure,  has  no  capacity  of 
imparting  kinetic  energy  to  the  exterior  flow. 


5.  Re-entrant  Jet  Model 

We  proceed  to  consider  the  effect  of  the  cavitation  number  a 
(defined  by(18))on  the  unsteady  flow  when  the  cavity  of  the  basic  flow  is 
finite  in  size.  An  additional  degreee  of  freedom  achieved  in  this  group  of 
flow  problems  is  that  the  cavity  volume  can  now  be  changed  arbitrary  by 
prescribing  a  flow  source  at  infinity.  In  order  that  the  basic  steady  flow 
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is  tractable  to  analysis,  resort  may  be  made  to  various  cavity-flow  models 
known  in  literature,  such  as  the  Riabouchinsky  model  [12],  the  re-entrant 
jet  model  [l],  the  wake  model  proposed  independently  by  Joukowsky  [13], 
Roshko  [14]  and  Eppler  [15],  and  the  modified  wake  model  introduced  by 
Wu  [16],  In  each  of  these  flow  models  an  artifice  of  some  sort  is  intro¬ 
duced  to  admit  o  as  a  free  parameter,  and  to  replace  the  actual  wake 
flow  of  a  real  fluid  by  a  simplified  model  within  the  framework  of  an 
equivalent  potential  flow.  It  has  been  found  that  in  a  positive  range  of  a 
up  to  moderate  values  of  order  unity,  the  agreement  between  these  models 
may  be  considered  very  close  (see  Wu  [5]).  Furthermore,  the  validity 
of  these  models  in  predicting  the  hydrodynamic  forces  acting  on  the  body 
has  been  supported  by  experimental  observations.  The  purpose  of  the 
following  several  sections  is  to  make  a  comparison  between  the  resultant 
unsteady  flows  when  different  models  are  used  for  the  basic  cavity  flow, 
for  such  a  task  should  be  of  fundamental  value  in  the  study  of  unsteady 
cavity  flows. 

For  simplicity  the  basic  flow  is  taken  to  be  that  past  a  flat  plate 
held  normal  to  the  stream  of  unperturbed  velocity  U  and  pressure  p^  , 
forming  a  finite  cavity  with  a  prescribed  cavity  pressure  Pj.(< 

According  to  the  re-entrant  jet  model,  the  free  streamlines  eventually 
reverse  their  direction  at  the  rear  part  of  the  cavity,  forming  a  re¬ 
entrant  jet  which  disappears  on  another  "Riemann  Sheet"  (see  figure  2). 
Due  to  the  assumed  symmetry  of  the  flow,  it  suffices  to  consider  only 
the  flow  in  the  left  half  physical  z-plane.  To  save  writing,  both  the  half 
plate  length,  i/2,  and  the  constant  speed  (see  (18)  )  on  the  cavity 
surface  will  be  normalized  to  unity. 

It  is  convenient  to  introduce  the  variable 

log  =  log  :=“  =  log  •^  +  i  e  (51) 

where  f^  is  the  complex  potential  of  the  basic  flow,  q^  is  the  flow 
speed,  and  B  the  flow  inclination  with  respect  to  the  x-axis.  The  part  of 
the  flow  under  consideration  in  the  z,f^,  auid  D-planes  is  shown  in  figure  2. 
By  applying  the  Schwarz-Christoffel  transformation,  the  flow  region  can  be 
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mapped  conformally  into  the  upper  half  (-plane,  with  the  point  of  infinity 
I,  front  stagnation  E,  plate  edge  A,  jet  infinity  B  and  the  rear  stag¬ 
nation  C  corresponding  to  (  =  «c,  0,  1,  b,  and  c  respectively.  The 
required  transformation  is  given  by 

3^  ^  ’  <52) 


dO 

3?  " 


_ B _ 

ta-c)[  ((-i)a-b)]^ 


(53) 


where  A,  B  are  two  coefficients.  From  the  local  behavior  of  12  at 
(  =  0  and  c,  we  find  the  relations 

B  =  -  ^  cb^  .  b  =  c(c-l)(c-i)’‘  •  (54) 

Integrating  (53),  and  by  making  use  of  (54),  we  obtain 


dz  Q 

W  =* 
o 


[(c-b)((-l)/f[c-l)((-b»l^-H 

l(c-b)(t.i)y8c.i)(c-b]ji*-i 


f[((-b)4b((-l)n^-H 


[ft-b)y!b((.i)|]KiJ 


(55) 

As  I  i  I  —  00  ,  dfjdz  — *  U  =  (l  +  <y)^.  From  this  condition  and  (55)  it 
results 


(2+a)/<T  =  ^  (c-1)^  c^*  (c-i)"*/*  .  (56) 


Finally,  by  integrating  (55)  and  (52)  to  obtain  z  =  z((),  the  coefficient  A 
is  determined  by  the  plate  length,  giving 

A  =(b-l)^*  /  [K(l)  -  K(o)]  .  (57a) 


where 
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K(e)  =  (2c.b)[b4(i-e)]^  +  (2c.i)[e(b.e)]^+(c-  \  b)fin"‘ 

+  b^[b(b-  I  )  +  c(2-b)]  sin"*(2$-l)  +  (b-l)^*(c-b)8in'‘ 


(57b) 


This  completes  the  zeroth  order  solution  for  prescribed  a.  In  particular, 
on  the  half  plate,  y=0,  -l<x  <0  (or  =  0+,  0  <  £  <  1),  we  have 


*(6)  =-  [K(e)  -K(o)l  /  [K(l)  -K(o)] 


(58a) 


\  ((c-b)(i.4)/i(c-i)(b-e)B^+i  1 

1  [(b-£)/!b(l-e)fl^+l  1 

J  i  ^ 

‘  [(c-b)(i.e)/«c.i)(b-eiw*-i  , 

t(b-e)^(i-©)]^-i 

This  expression  is  needed  for  the  first  order  solution. 

When  a  sudden  acceleration  of  magnitude  b^  is  applied  normal  to 
the  plate,  the  boundary  conditions  corresponding  to  (32)  can  be  written 


g'*’  +  o' 

=-2i  b^  dx/dg 

for 

0  <  e  <  1  . 

(59a) 

G'*’  -  G" 

=  0 

for 

1  <  e  <  b  , 

(59b) 

g'*'  +  g’ 

=  0 

for 

.«  <  1  <  0 

and  b  <  5  ^  • 

(59c) 


The  last  condition  (59c)  expresses  the  assumption  that  the  perturbed 
flow  preserves  the  basic  flow  symmetry  about  the  y-axis,  that  is, 

V  =0  on  11  =  0,  for  ^  <  0  and  £  >  b.  The  above  conditions  are  ex¬ 
pressed  in  the  form  different  from  (32).  The  solution,  however  can  be 
written  down  directly  by  applying  the  same  principle.  It  is  obvious  that 
H(i)  =  [  (C-l)((i-b)  ]  *,  defined  with  branch  cuts  from  (i  =  -oo  to  1  and 
from  b  to  00  so  that  H  — *■  (as  |  ( |  oo  ,  0  <  arg  I  <  ir,  is  a  homo 
geneous  solution  of  the  present  problem.  Therefore  the  general  solution 
can  be  written 
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n 


-00 


(60) 


where  dx/d|  is  given  by  (58b)  and  are  arbitrary  real  constants.  To 

determine  c^,  we  note  first  that  the  pressure  must  be  finite  at  {^  =  0, 

hence  =  0  for  n<  0.  Furthermore,  we  note  that  the  first  term  in  (60) 

is  of  order  of  4'*  as  |C|“*  «  .  which  implies  that  the  behavior  of 

df^dz  at  large  distances  is  determined  by  the  last  term  in  (60).  By 

applying  condition  (30)  and  using  the  symmetry  property  of  the  flow,  we 

obtain  c^  =  Q  /2'tr,  c  -  c  =  .  .  .  =  0,  where  Q  is  the  source  strength 
O  1  1  2  1 

defined  in  (24). 


The  integral  in  (60)  can  be  integrated  explicitly.  For  the  deter¬ 
mination  of  hydrodynamic  forces,  however,  only  the  real  part  of  G  on 
the  plate  (where  r]  =  0,  0  <  4  <  1)  is  needed.  The  final  result  is 


R,o  =  .  i 

Ik(1).K(o)1 


A  +B  (1-1) 

i  I 


tjb  (t-l)[  (2c-b-2)({-l)-c(b.l)l  (b-l){t(b:{y2j(b.l){(b-U 


[e(b-4)i= 


i  [(l-4)(b-|)I’^ 


(61) 


where 


=ir|^b(c-l)+  ^1  +  '^  (b*-2b+2c)  (-5  +sin'‘  j  -  Jb(b-l)(b-2c)  , 

B^  =/^(b-2c)^-2  +  sin  ‘  -ir(2c-l)  . 

By  using  (40b)  and  the  symmetry  property  of  the  flow,  the  first  order 
normal  force  is  given  by 

N  =-2p  (  x(e)  {RiGldg  ,  {62a) 

^  ‘'o 


where  x(g )  and  Ri  G  are  given  by  (58a)  and  (61).  The  result  can  be 
expressed  in  terms  of  nondimensional  parameters  as 
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in  which  the  coefficients  and  Fq  are  functions  of  the  cavitation 

number  a.  Analytic  evaluation  of  the  above  integral  is  too  tedious  to  be 
practical.  The  numerical  computation  of  these  coefficients  has  been 
carried  out  with  an  IBM  7090,  the  final  result  is  plotted  versus  a  in 
figures  6  and  7  to  show  the  effect  of  (i)  the  finite  cavity  size  and  (ii)  the 
change  in  the  cavity  volume.  From  this  result  several  salient  features  of 
significance  may  be  pointed  out. 

First,  it  is  noted  that  as  a  — ►  0,  the  value  of  F^  tends  to 

0.8448,  which  is  the  limit  of  Kirchhoff-Helmholtz  case  (see  equation  (42) ). 

From  small  to  moderately  large  values  of  o,  F^  increases  very  slowly 

1 

with  increasing  a  as  compared  with  the  rate  of  increase  of  the  zeroth 
order  drag,  which  increases  approximately  with  the  factor  (1+cr)  (see, 
e.g. ,  Gilbarg  [l]  ). 

Another  point  of  interest  is  that  Fq  0  rather  rapidly  as 

9 — *  0,  this  limit  being  independent  of  so  long  as  is  finite. 

This  result  shows  that  the  effect  on  the  drag  force  of  removing  fluid  at 

infinity  is  insignificant  when  the  cavity  is  sufficiently  long.  Furthermore, 

it  shows  that  the  limit  of  the  solution  as  9  — -  0  is  non-uniform  with 

respect  to  Q  since  the  solution  in  the  Kirchhoff  case  does  not  exist 
1 

unless  Q  =  0. 

1 

6.  Riabouchinsky's  Model 

The  essential  feature  of  the  Riabouchinsky  model  is  the  introduct¬ 
ion  of  an  appropriate  image  body  downstream  of  the  real  body  so  that  the 
free  boundaries  of  the  cavity  are  connectedby  this  pair  of  solid  boundaries. 
Let  us  apply  this  model  to  consider  the  cavity  flow  past  a  flat  plate  set 
normal  to  the  stream,  the  flow  in  the  physical  z-plane  is  shown  in 
figure  3.  Again,  as  in  the  problem  stated  in  the  previous  section,  the 
unperturbed  velocity  and  pressure  are  U  and  respectively,  the 
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cavity  pressure  is  p^,  corresponding  to  the  cavitation  number  a.  Also 
the  half  plate  length,  i/2,  and  the  constant  speed  on  the  cavity  will 
be  normalized  to  unity.  Furthermore,  due  to  the  assumed  symmetry, 

only  the  left  half  z-plane  of  the  flow  need  be  considered. 

For  the  present  case  it  is  convenient  to  denote  by  C  the  complex 
velocity  potential.  We  further  introduce  an  auxili^i^y  complex  variable 
T  defined  by  . 

where  w^  =  d^/dz  is  the  hodograph  plane  of  the  basic  flow.  The  flow 
field  under  consideration  in  the  z-plane  and  conformal  mapping  planes 
w^,  T,  are  shown  in  figure  3.  By  the  assumed  synunetry  of  the  flow, 
we  may  choose  the  potential  at  the  front  and  rear  stagnation  points  to  be 
^  =  -n  and  n  respectively,  and  =  -m,  4^  =  m.  We  further  define 
k'  and  k  by 

•  k  =(l-k'*)^  =  (2<T  +  a*)^/(l  +a)  .  (64) 

Then,  at  z  =  oc  ,  we  have 

Wj= -ijk'  .  Tj  =  ^  {k''^-k'^}  (65) 

The  upper  half  r-plane  is  mapped  into  the  upper  half  C-piS'ne  by  the 
Schwarz-Christoffel  transformation 

&  =  2njir’  T[4k'T*  +  (l-k')*]"^  .  (66) 


From  the  local  conformal  behavior  of  ^T)  at  the  point  D,  i.e. 

=  m,  Tq  =  1,  we  find  the  relationship 

m  =  2njir/(l  +  k')  .  (67) 


Let  us  introduce  another  auxiliary  variable  u  =  fi  +  iy? defined  by 

(68) 


w 


=  =  -i  ^ 

3z  dn  u 


where  dnu  is  the  Jacobian  elliptic  function,  delta  amplitude  of  (o,  with 
modulus  k.  In  the  following  analysis  the  conventional  notations  (see  [17]) 
for  the  elliptic  functions  and  integrals  will  be  used  without  specification  and  the 
modulus  k  will  always  be  omitted  to  save  writing.  By  substituting  (68) 
euid  (63)  into  (66)  it  gives 
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C  =  i 


n  dn*  fa)-k' 
“K"  anu  cno) 


(69) 


From  (68)  and  (69)  we  deduce  that  on  both  the  front  and  the  image  half 
plates  -  1  <  X  <  0, 


*=  .  0<t‘<7  . 

(70) 


Evaluating  the  above  result  At  fj.  -  K/2,  which  corresponds  to  x  =  -1, 
we  obtain 


n 


2 


E  -  k'K  + 


d-k')* 

T+R'" 


(71) 


which  completes  the  necessary  calculation  for  the  basic  flow. 

Due  to  the  presence  of  the  image  plate,  an  additional  assumption 
is  needed  for  this  model  in  the  study  of  its  small  time  behavior.  It  is 
given  that  at  t  si  0  a  sudden  acceleration  b^  is  applied  on  the  front 
plate  directed  in  the  positive  y-direction.  For  small  t  >  0  the  image 
plate  may  be  assumed  to  move  in  the  y-direction  with  speed 

V  = -pb^  t  +  0(t*)  ,  (72) 

where  p  is  an  unknown  constant.  To  determine  this  unknown  constant 
P,  we  shall  assume  that  on  the  image  plate  the  jump  in  the  drag  due  to 
the  suddenly  applied  acceleration  is  zero.  The  physical  significance  of 
this  assumption  may  be  explained  as  follows.  It  has  been  pointed  out  by 
Wu  [5]  that  the  image  plate  in  the  basic  flow  may  be  regarded  as  a 
means  to  represent  the  energy  dissipation  in  the  wake  flow  of  a  real 
fluid.  In  fact,  in  a  frame  of  reference  at  rest  with  respect  to  the  fluid  at 
infinity,  the  work  done  by  the  moving  image  plate  is  negative  and 
numerically  equal  to  the  work  done  by  the  real  plate  since  the  total  force 
on  the  pair  of  plates  vanishes.  This  negative  work  done  by  the  image 
plate  therefore  corresponds  to  the  mechanical  energy  removed  from  the 
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system  in  unit  time  as  there  is  no  other  means  of  dissipating  energy  in 
potential  flow.  Now  in  the  unsteady  motion,  it  is  conceivable  that  the 
rate  of  dissipation  in  the  wake  (of  the  actual  flow)  cannot  be  affected  at 
small  time.  This  implies  that  the  initial  change  of  momentum  at  the 
image  plate  must  vanish.  It  is  this  physical  meaning  that  underlies  the 
above  assumption.  For  small  t  >  0,  the  boundary  condition  on  the  part 
of  the  %  -axis  (ti  =  0^)  corresponding  to  the  image  plate  can  be  written 
by  using  (72)  as 

ImG=  pbR<(^j, 

where  G  =  df^/d4.  Or,  by  applying  the  continuation  (31),  we  have  on  the 
part  of  the  4  -axis  corresponding  to  the  image  plate  the  condition 

G‘‘'  +  G‘=  Zipb^Rf^^j.  (73) 

Referring  to  conditions  (32)  and  (73)  and  using  (68),  the  boundary 
conditions  of  this  problem  are 

onG-  =  .ab_R/(^). 

G"*"  -  G"  =  0  for  -m  <  §  <  m 

G+tG-=  afbRl(ll).  a^^^£orm<e<n 

G^  +  G'  =  0  for  -oo<4<-n  ,  n<|<oc  ; 


where  m  and  n  are  given  by  (67)  and  (71). 

Consider  the  function  h(4)  -  m*  with  the  branch  cuts  from 

-oc  to  -m  auid  from  m  to  oo  ,  and  hd) — ►  ^  as  |l| — ►  oo  , 

0  <  arg  ^  This  h(^)  is  obviously  a  solution  of  the  corresponding 

homogeneous  Hilbert  problem  and  satisfies  (31),  Therefore,  by  using 
(33),  the  required  solution  may  be  written 


G  (;)  =  - 


b  k 
1 


1 

i-I 


fZ 

•j;  -m 


n  -m 


e*-m 


i  “-  I 


de 


-4” 


2  cn^ 

IS 


dt 


+ 


P(0 

C*-m* 


(74) 
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ec 

Where  P(C)  =  ^  and  are  arbitrary  real  constants.  Again, 

-  oc 

since  the  pressure  at  4=0  should  be  finite  and  the  first  term  in  (74) 
is  of  order  of  4  *  as  |4|  ^  ,  we  have,  by  applying  condition  (30) 

and  the  symmetric  property  of  the  flow,  P(4)  =  Q  /2it.  By  using  (69), 
the  variable  4  i^i  (74)  can  be  transformed  into  the  variable  fi,  giving 


G(4)  = 


b^(l-k')n 
rrjk'  J4*.m* 


(75) 


After  integrating  the  above  expression  and  by  using  (69)  we  deduce 

that 


where  the  upper  sign  is  used  for  the  front  plate,  the  lower  sign  for  the 
image  plate,  w*  denotes  the  function  which  corresponds  to  the  Cauchy 
principal  value  of  the  integral  representation  of  the  function  w  (for  the 
definition,  see  [17]),  and 


28 


A 

z 


1-k' 

T+lc' 


E+k'K+d+k') 

IT  “k')^ 


B 

2 


^  j.  «•  -1  1-k' 

■z  r+i?  +  t+f 


{76b) 


{76c) 


By  applying  (40b),  (70),  (76a)  and  (66),  we  obtain  the  normal  force  N  on 

1 

the  front  plate  as 

N  =  2p  {bJ(l+P)I^  +  (l-p)!^]  +  0^13}  .  (77) 

and  the  normal  force  (in  the  positive  y-direction)  on  the  image  plate 
as 


Dy  =  2p  {bj(l+p)l^  -  (l.p)lJ+Q^I^  } 


(78) 


where 


I  -  x(u)  {1+k')  JiT  *  .  K  (dnV-k') 

‘  irXF  Jo  ^  n  2  7  ^2^ 


+  (dn*pt-k' 


'  sin~^  -m}i  ■ 

.  '  ) 

.  snV 

^  Jl+k' 

sn*/i  ' 

.i3 


.iSliL  ,  /.In-I  ‘  ,  52i\] 

bicn^fl  \  jl+k'  cn*p  /J 


dn*/icn* 

^  r 


dM  > 


{79a) 


B 


dn' 


Hii.')  sn,.cn^  +  ^ 


dnV  y 


(l-k')dnV 


.  log  jEi 

cnfi~  Jk'  anfi 


dfi 


{79b) 


^3  =  “7¥ 


‘-'‘1  x(rt 


iifi 


dfi 


{79c) 
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and  x(^)  is  given  by  (70). 

Now  we  apply  the  condition  that  on  the  image  plate  the  change  in 
the  drag  due  to  the  applied  acceleration  is  zero,  giving 

2pb^[  (l+p)I^  -  (l-P)lJ  =0  , 
or 

P  =-r^  • 

1  2 

Substituting  (80)  into  (77),  we  obtain  the  drag  on  the  front  plate  as 


N  = 

t 


(81) 


This  result  can  also  be  expressed  in  the  following  non-dimensional  form: 


where 


(82a) 


Fq  =  -21,  .  (82b) 


The  integrals  and  are  computed  numerically  with  an  IBM  7090 

amd  the  final  result  of  F,  and  F.  are  shown  in  figures  6  and  7  as  a 

comparison  with  the  results  of  other  flow  models.  It  is  further  noted  that 

as  <T  — ►  0,  F,  reduces  to  (42)  of  the  Kirchhoff  case  and  F,..  — ►  0. 

D  U 

1  1 

7.  A  Wake  Model  for  an  Oblique  Plate  with  a  Finite  Cavity 

Thus  far  we  have  considered  the  accelerating  motion  of  an  inclined 
plate  in  Kirchhoff  flow,  and  of  the  finite-cavity  flow  past  a  plate  broad¬ 
wise  to  the  stream.  Application  of  either  the  re-entrant  jet  model  or 
the  Riabouchinsky  model  to  an  oblique  plate  with  a  finite  cavity  formation 
leads  to  very  complicated  analyses.  The  task  is  considerably  simplified, 
however,  if  we  adopt  a  modified  wake  model  recently  proposed  by  Wu[  16] 
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to  describe  the  basic  steady  flow.  The  purpose  of  this  section  is  to 
determine  the  effect  of  acceleration  on  the  hydrodynamic  forces  with  the 
cavitation  number  a  and  the  angle  of  attack  a  as  two  free  parameters. 

The  basic  flow  is  taken  to  be  a  uniform  stream  of  infinite  extent 
impinging  on  a  flat  plate  at  an  incidence  angle  a,  to  which  a  finite  cavity 
is  attached.  According  to  this  modified  wake  model,  the  incoming 
stagnation  streamline  branches  off  the  plate  at  the  leading  edge  A  and 
the  trailing  edge  B,  forming  two  free  streamlines  ACI  and  BC'I  which 
are  assumed  to  become  asymptotically  parallel  to  the  main  flow  at  down¬ 
stream  infinity  (see  figure  4).  The  pressure  on  the  parts  AC  and  BC 
of  the  free  streamlines  is  assumed  to  take  the  constant  cavity  pressure 
p^,  and  the  space  within  the  closed  curve  ACC'BA  is  regarded  to 
represent  approximately  the  cavity.  The  space  in  between  the  free 
streamlines  Cl  and  C'l  represents  a  crude  model  of  the  dissipating 
wake,  along  its  boundary  the  uniform  stream  conditions  are  eventually 
restored  at  downstream  infinity.  The  flow  outside  this  infinite  wake  strip 
is  assumed  to  be  irrotational.  The  locations  of  the  points  C  and  C  are 
determined  with  two  assumptions;  the  first  is  that  both  the  velocity 
potential  and  the  flow  inclination  at  C  and  C  are  equal,  and  the  second 
assumption  is  the  so-called  "hodograph-slit  condition"  that  the  free 
streamlines  Cl  and  C'l  form  a  slit  of  undetermined  shape  in  the  hodo- 
graph  plane.  With  these  two.  additional  assumptions  the  whole  flow  field 
outside  the  wake  is  then  completely  determined.  For  the  convenience  of 
subsequent  application,  the  solution  of  the  basic  flow  is  reproduced  briefly 
in  the  following.  The  plate  length  t  and  the  constant  speed  q^  on  AC 
and  BC  are  again  normalized  to  unity. 

The  flow  in  the  physical  z-plane,  the  complex  potential  f^-plane, 
and  the  hodograph  \^-plane  are  shown  in  figure  4.  The  subscript  of  w^ 
will  be  omitted  for  simplicity.  We  further  introduce  the  parametric 
(-plane  defined  by 


(  =  4  (w  +  w‘*  ) 


(83a) 
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or 

w  =  t  -  (I*-  1)^  ,  (83b) 

in  which  the  function  1)‘  is  defined  with  a  branch  cut  made  between 

the  points  4  =  -1  and  1  so  that  (4*  -  1)^  “*  4  as  |4l'~**o-  At  the 
point  of  infinity,  the  complex  velocity  taJces  the  value 

.  i 

w=W=Ue*^®  ,  U=(l+ff)'“  .  (84) 


The  corresponding  value  of  4  is 

^00=  2(U'‘e'"  +  Ue‘"“)  .  (85) 

Since  Im  =  0  on  the  entire  real  4-sxis,  the  complex  potential  fJ4) 
can  be  continued  analytically  into  the  lower-half  4'Plsne  by 

£^(7)  =  •  (86) 


Now  from  the  asymptotic  behavior  of  the  streaumlines  =  const,  near 

the  point  4  -  4gg  >  it  is  evident  that  f^  must  have  there  a  simple  pole. 

Furthermore,  from  the  local  conformal  behavior  of  f  at  f  =  0,  it  is 

o  o 

obvious  that  f^  =  0(4~*  )  as  |  4  I  *  •  Therefore  the  solution  must  be 
of  the  form 


f 

o 

where  A  is  a  real  constant. 


9  i -  (87a) 

Or,  expressing  in  terms  of  w  by  (83), 


_ A  w^ _ 

o~  (w-W)(w-'W')(w-W*)(w-V*) 


(87b) 


The  z-plane  is  determined  by  integration  of  dz/df^  =  1/w,  giving 


fo(w) 


w 


-  ( 


w 


+  i  B  (  ^  -  )  1^^  log  (w-W)  -  W  log  (w  -  ^  )j 

^  log  (w  -  W  )  -  W  log  (w  -  ^  )  j 


(88a) 


z  +  a  = 
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where  the  constant  B  is  related  to  A  by 


g'  =  2  (  ^  -  U)  sin  a 


g  +  U)*  -  (2  cos  ar)*  j  , 
Finally,  the  constant  A  is  determined  by  the  plate  length  as 


(88b) 


{89a) 


A  =  [^(  g  +  U)*  -(2  cos  a)^j  /k  . 

,  (U‘'+U)^+(2co8o)2  .  ir(U‘*+U)  .  (U‘'+U)"-(2  cos  a)*  ') 


(89b) 


For  the  unsteady  motion  we  shall  confine  ourselves  to  the  simple 
case  of  constant  acceleration  so  that  for  any  point  z  of  the  plate, 

dz/dt  =  (aj+  ib^)t  .  (90) 

Then  the  boundary  conditions  of  this  problem  (see  equation  32)  become 
G'*'  +  G"  =  -  2i  b  R1  (dz/d&)  for  |g  1  >  1  , 

+  -  *  i  I 

G^  -  G  =  0  for  ||  I  <  1  , 

where  G  =  df  /d^.  In  the  first  condition,  dz/d^  =  w  *  df^/dC  can  be 
deduced  from  (83)  and  (87). 

Similar  to  the  general  solution  of  the  problem  with  the  boundary 
conditions  given  in  (32),  the  general  solution  of  the  above  boundary  value 
problem,  stated  in  (91),  can  be  written 


G(&)  = 


ir(;*-l)^ 


di  +(t*-l)"^P(t) 

(92) 


where  the  function  (^^  -1)^  is  defined  in  the  entire  (-plane  with  branch 

cuts  from  -oo  to  -1  and  from  1  to  oo  along  the  |-axis  so  that 

((*  -  1)^—*  (as  1(1  — *  00  ,  for  0  <  arg  (  <  xr  (note  that  the  branch  cut 
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for  the  function  (4^  '  1)^  for  this  perturbation  part  is  different  from  that 
for  the  basic  problem  as  described  after  (83b)).  The  arbitrary  function 
P({^)  is  real  on  the  real  4-&xis,  and  hence  can  be  expanded  into  the 
Laurent's  series 

oc 

n=  -oc 

Since  z  'as  |z  |  — ►  oc,  it  is  necessary  to  have  c^  =  0  for 

n  <  -  2  in  order  that  the  perturbation  velocity  w  vanishes  at  infinity. 
Furthermore,  we  must  impose  c^^  =  0  for  n  >0  if  we  require  the 
pressure  at  1^1  =oo  to  be  finite.  Finally,  application  of  (30)  yields 


P(&)  =  +c_^(4-;^)''  .  (93a) 

where  use  has  been  made  of  (84)  and  (85). 

Transforming  the  variable  ^  in  the  above  solution  into  w  by  (83) 
we  obtain 


2b  w* 

I 

ir(l-w*) 


Q  w^ 

1 

■ir(l-w*) 


1-W*  j.  l-W'^ 

(w-W){wW-l)  (w-'W^(wW-l) 


(94a) 


where 

/  X  1  1-u^ 

=  u  -air  ({r-w)rdw~r)  • 

In  particular,  on  the  rigid  plate,  -l<w<l,  RfG  =  8^^/8S,  and 


(94b) 


2w^ 

it(1-w*) 


b^(P.V.) 


r  1  Q 

J  ^  g(w,u)du - ^ 


1-W* 

(w-W)(wW-l) 


+ 


1-W^ 

(w-W)(wW-l) 

4 


(95) 
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where  (P.  V. )  denotes  the  Cauchy  principal  value  of  the  integral. 

The  normal  force  N  acting  on  the  plate  due  to  the  acceleration 
is  given  by 


N 

1 


=  P 


x(w) 


(96) 


where  8^  /8|  is  given  by  (95),  and  d^/dw  can  be  obtained  from  (83). 
As  in  the  previous  cases,  we  write 


then,  from  (96), 

Tb  =  I  (P.V.)  x{w)dw 


(97) 


(98) 


=  .  i  r 

\  ^  J-l 


1  -W* 

■(\^TW)Xw'wnT 

k 


i-w* 


(w.W)(wW-l) 


dw  . 


(99) 


These  integrals  have  been  computed  numerically  with  an  IBM  7090,  the 
final  result  of  and  Fq  is  shown  in  figures  8  and. 9  for 

o  o  ‘  o  * 

a  =  90”,  75”  and  45”.  The  result  of  the  special  case  a  =  n/ 2  is  also 
compared  with  the  other  flow  models  in  figures  6  and  7.  For  the  special  case 
a  =  Tt/Z,  we  deduce  from  (98)  the  following  expression  for  small  a 

F^^  =  0.8448  [  1  +  0.067 a^l-o)  +  (a'^.)] 

1 

which  is  the  special  case  already  noted  by  Wu  [6]. 


(100) 
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8.  Conclutlon 

In  the  foregoing  sections  we  have  treated  the  accelerating  motion 
of  an  inclined  plate  in  Kirchhoff  flow  and  in  the  finite-cavity  flow,  and 
that  of  a  pure-drag  plate  in  Riabouchinsky  flow  and  in  the  re-entrant  jet 
flow.  From  these  results  we  make  conclusions  by  pointing  out  the  follow¬ 
ing  features  of  interest. 

(i)  Effect  of  the  cavitation  number  a -  For  0  <a  <  1,  the 

force  coefficient  due  to  the  normal  acceleration  evaluated  by  using 

1 

three  different  cavity-flow  models  is  nearly  identical,  tending  to  0.8448, 

the  Kirchhoff-Helmholtz  limit,  as  a — ►  0,  and  being  very  insensitive  to 

a  in  this  range.  As  a  increases  further  from  1,  increases  very 

1 

slowly  compared  with  the  rate  of  increase  of  the  steady  state  drag  co¬ 
efficient,  which  has  the  factor  (1  +  a)  approximately.  For  a  fixed  o 

in  this  region,  the  value  of  by  the  wake  model  is  higher  than  that  of 

1 

the  Riabouchinsky  model,  which  is  in  turn  higher  than  that  of  the  re¬ 
entrant  jet  model.  This  result  indicates  that  in  all  practical  ranges  of  a, 
it  is  unimportant  to  consider  the  effect  of  o  for  the  small-time  limit  or 
the  large  frequency  limit. 

(ii)  Effect  of  incidence  a -  For  moderate  and  large  o,  the 

force  coefficient  F^  increases  faster  with  increasing  o,  the  smaller 

1 

is  the  value  a.  This  result  should  be  expected  on  the  physical  standpoint 
since  there  contains  more  liquid  in  the  circular  cylinder  spanned  by  the 
plate  at  smaller  a.  This  effect  of  a,  however,  becomes  insignificant 
for  a  less  than  order  unity. 

(iii)  Effect  of  change  in  cavity  volume -  The  force  coefficient 

Fq  due  to  the  change  in  cavity  volume  displays  wider  differences  be- 

tween  different  flow-models.  All  these  models  show,  however,  that 
Fq  — *  0  rather  rapidly  as  a — ►  0,  saying  that  the  effect  on  the  drag 

force  of  removing  fluid  at  infinity  is  insignificant  when  the  cavity  is 
sufficiently  long.  This  analysis  also  justifies  the  assumption  that  p  is 
bounded  at  infinity  introduced  in  some  of  the  previous  works  [4,6,7,8,9,18] 
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in  the  case  of  unsteady  flows  with  infinite  cavities.  As  shown  by  this 
analysis,  this  assumption  actually  turns  out  to  be  a  necessary  condition 
for  the  existence  of  the  solution  in  the  case  of  a  =  0. 
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Figure  5.  Normal  force  coefficients  for  inclined 
lamina  in  Kirchhoff  flow. 


Figure  6.  Normal  force  coefficients  for  the  sudden  acceleration 
of  a  flat  plate  placed  normally  to  the  incoming  stream. 


Figure  7.  Normal  force  coefficients  due  to  a  source  at  the  point  of  infinity 
for  a  flat  plate  placed  normally  to  the  incoming  stream. 
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Figure  9.  Normal  force  coefficient  due  to  a  source  at  the  point  of  infinity 
for  a  flat  plate,  at  an  angle  of  attack  a,  in  the  wake  model. 
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